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We perform a first global exploration of the Constrained Next-to-Minimal Supersymmetric Stan- 
dard Model using Bayesian statistics. We derive several global features of the model and find that, 
in some contrast to initial expectations, they closely resemble the Constrained MSSM. This remains 
true even away from the decoupling limit which is nevertheless strongly preferred. We present ensu- 
ing implications for several key observables, including collider signatures and predictions for direct 
detection of dark matter. 



I. INTRODUCTION 

Effective low-energy supersymmetry (SUSY) has many attractive features and is widely expected to provide a more 
complete description of phenomena at and above the electroweak scale than the Standard Model (SM) of electroweak 
and strong interactions. In many SUSY models, gauge coupling unification can easily be achieved, unlike in the SM 
or non-SUSY versions of those models. Moreover, SUSY, when supplemented by R-parity (or matter parity), offers 
an attractive candidate for resolving the dark matter (DM) problem in the Universe. 

On the other hand, SUSY, being a global symmetry, allows for a whole multitude of possible effective models which 
otherwise would suffer from the well-known hierarchy and fine-tuning problems. In the past, most phenomenological 
studies focused on the Minimal Supersymmetric Standard Model (MSSM) - a supersymmetrized version of the SM [l[ . 
More recently, a constrained version of the MSSM (CMSSM) 0, which includes a minimal supergravity (mSUGRA) 
model, has become more popular by virtue of its relative simplicity and a small number of free parameters. This is 
achieved by relating at the unification scale the soft masses of the MSSM gauginos to a common value mi/2, those of 
the scalar partners of SM fermions to mo, the tri-linear terms to Aq, in addition to tan/3 - the ratio of v.e.v.'s of the 
neutral components of the two Higgs doublets. 

One puzzling and unsatisfactory feature of the MSSM is the so-called /i-problem [3|: the Higgs/higgsino mass 
parameter is SUSY-preserving but, on phenomenological grounds, it is expected to be of the same order as soft SUSY 
breaking masses, [i ~ A/susy — 1 TeV. Various solutions have been suggested, for example 

A model that solves the /i-problem of the MSSM in a simple way is the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM) Q. In the NMSSM one adds a singlet chiral superfield S. The explicit fi term is absent and the 
superpotential has only dimensionless parameters and therefore the only new effective scale is the one of the soft 
breaking terms AfsusY- The fi parameter is generated dynamically through the v.e.v. of the spin-0 component of the 
singlet superfield. 

At the phenomenological level, the presence of additional fields, namely an extra CP-even and CP-odd neutral Higgs 
bosons, as well as a singlino component of a neutralino, leads in general to a richer and more compl ex p henomenology p, 
0) HI, H, EH , as well as cosmology, in particular with respect to the domain wall problem [ll|, Q2, [H| . 

In analogy with the CMSSM, successful gauge coupling unification in SUSY has provided motivation for considering 
a constrained version of the NMSSM (CNMSSM) [IJ, fl5fl . which we will define below. Extensive phenomenological 
investigations of the CNMSSM have been carried out in j6[. 

On the other hand, the enlarged set of parameters makes a full exploration of the CNMSSM even more challenging 
than in the case of the CMSSM. Traditional techniques of sampling slices of the parameter space provide limited 
information and are inadequate in a number of other aspects, for example in fixing relevant SM parameters which 
may have much impact on the outcome, as shown in (lfjl. Il7l. Il8j. 

More recently, it has been demonstrated that a Markov Chain Monte Carlo (MCMC), or some other, scanning 
technique, coupled with Bayesian statistics, can very efficiently pro be the whole parameter space and thus allow one 
to derive global properties of the model under investigation [lg. Tl7l [H, [l9| . 

Over the last few years, several studies using the Bayesian approach have been performed of the CMSSM [l6|, 03, 
[IIGjI, the Non-Universal Higgs Mass Model (NUHM) 0, the MSSM [HSU and large volume string [H|. 

One advantage of the Bayesian approach is that it allows quantitative model comparison, using so called Bayes 
factors to see which model best fits the data, be it selecting the sign of n in the CMSSM to picking a class of SUSY 
breaking. 
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In this paper we perform a Bayesian analysis of the CNMSSM. We explore very broad ranges of the CNMSSM 
parameter and apply all most important experimental and cosmological constraints, including all collider limits, the 
branching ratio of b — > 57, the difference 5(g — 2) M between the experimental and SM values of the magnetic moment 
of the muon, the LEP limits on sparticle and Higgs masses and the 5 year WMAP limit on the relic abundance fl x h 2 
of the lightest neutralino assumed to be the dark matter. A full list of constraints used and the exact numbers used 
in the analysis will be given below. 

Our main finding is that, somewhat contrary to initial expectations, from the statistical point of view, most 
phenomenological and dark matter features of the CNMSSM of crucial importance for experimental tests closely 
resemble those of the CMSSM. In particular, singlino-dominated LSP only appears in a very limited number of cases 
that are not yet excluded by experimental bounds on the parameter space. Clearly, this will make it very challenging, 
although not impossible, to distinguish the models at the LHC and in DM searches, as we discuss below. 

The paper is organized as follows: in Sec. |TT] we define and overview the NMSSM and the CNMSSM. In Sec. IIIII 
we describe our the statistical approach. The results are presented in Sec. IIVI and we finish with our conclusions in 
Sec. El 



II. THE NMSSM AND THE CNMSSM 



The NMSSM superpotential contains a new superfield S which is a singlet under the SM gauge group SU(3) C x 
5?7(2)l x U(1)y- (We use the same notation for superfields and their respective spin-0 component fields for simplicity.) 

W = ey (Y u W u Q l u + Y d H l d Q j d + Y e H d V e) - ey A S H d W u + ^kS 3 , (1) 

where if J = (ff9, H d ), ffj = (ff+, ff°), i,j are SU(2) indices with £12 = 1, while A and k are dimensionless couplings 
in the enlarged Higgs sector. 

The superpotential in Eq. ([T]) is scale invariant, and the EW scale will only appear through the soft SUSY breaking 
terms in £ so ft , which in our conventions is given by 

— £ S oft = ttVq Q* Q + to 2 -, it* u + m 2 ^ d* d + m 2 ^ L* L + m 2 -, e* e 
+ m 2 Hd H* d H d + m% n H*H U + m% S* S 

+ ^ (A u Y u W u Q l u + A d Y d H d & d + A e Y e H d LP e + H.c.) 
-djXAxSHiHi + ^kA k S 3 + H.c. 

- - [M a A 3 A 3 + M 2 A 2 A 2 + Mi Ai Ai + H.c.) . (2) 

When the scalar component of S acquires a VEV, s = (S), an effective interaction [iH d H u is generated, with fi = As. 
In addition to terms from £ S oft, the tree-level scalar Higgs potential receives the usual D and F term contributions: 

Vd = - \H U \ 2 ) 2 + I \H\H U \ 2 , 

V F = |A| 2 {\H d \ 2 \S\ 2 + \H U \ 2 \S\ 2 + \e tJ H d Hi\ 2 ) + | K | 2 |S| 4 

-(e ij \K*H d HiS* 2 + K.c.) . (3) 

Using the minimization equations we can re-express the soft breaking Higgs masses in terms of A, k, A\, A k , 
l 'd = (H d ), v u = (ff°) (with tan/3 = v u /v d ), and s: 

m 2 Hd = - A 2 (s 2 + v 2 sin 2 0) - \ M z cos 2 P + As tan P ( KS + A \) > ( 4 ) 
m% u = - A 2 (s 2 + v 2 cos 2 0) + -M| cos 2(3 + As cot (3 [ks + A\) , (5) 

ml = - X 2 v 2 - 2n 2 s 2 + Xkv 2 sin 2/3 H — sin 2(3 - nA K s , (6) 

2s 

The boundary conditions at the grand unification scale Af^uT — 2 x 10 16 GeV are analogous to those of the 
CMSSM, with the exception of mj. n, M s and s are fixed by the minimization equations (dJ)-© which leads to five 
continuous free parameters of the CNMSSM: TO1/2, mo, Aq, tan/3 and A, in addition to sgn(/i). 
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The feature of not unifying m s with all the other soft scalar masses at mo gives one the necessary freedom to 
obtain, in the limit A — > 0, with As fixed, the CMSSM plus a singlet and a singlino field that both decouple from the 
spectrum, as discussed in [lfij]. For cosmological analyzes those extra particles can still play an important role but 
from the particle phenomenology point of view the model becomes indistinguishable from the CMSSM. Also in this 
limit a singlino LSP is excluded since, being completed decoupled, it can not annihilate into SM particles. 

We also present the neutralino sector since the lightest neutralino will, by assumption, play the role of dark matter. 
The mass term in the Lagrangian is given by 



^ SS = 4(*TAV*° + H.c., (7) 



with M.^,o given by a 5 x 5 matrix, 



M x o 



Mi -M z sin0 w cos/3 M z sin 6 W sin (3 

M 2 M z cos 9 W cos j3 -M z cos W sin/3 

-Mz sin 9w cos /3 Mz cos 9w cos /3 —As — Xv u , (8) 

Mz sin dyv sin /3 —Mz cos 8yy sin /3 — As — \v<i 

V —Xvu ~Xvd 2ks ) 

where M x (M 2 ) denotes soft the mass of the bino (wino) and 8w denotes the weak mixing angle. 

III. OUTLINE OF THE METHOD 

Following the discussion of the previous Section, in the CNMSSM the free parameters are given by 

9 = (m 1/2 ,m , A ,tanf3, A), (9) 

while we fix sgn(/x)= +1, which implies s > 0. Furthermore, without loss of generality we choose A > 28]. 

As the values of relevant SM parameters, when varied over their experimental constraints, have an impact on the 
observable quantities, fixing them would lead to inaccurate results. Instead, here we incorporate them explicitly as 
free parameters (which are then constrained using their measured values), which we call nuisance parameters ip, where 

V> = {Mum^m^ ; Qs H M ). (10) 

In Eq. (Tl0|) M t denotes the pole top quark mass, while the other two parameters: mb(mb) MS - the bottom quark mass 
evaluated at TOf, and a s (Mz) M - the strong coupling constant evaluated at the Z pole mass mz - are all computed 
in the MS scheme. Note that, in contrast to recent analyzes of the CMSSM [l6|, [13, Ha ], we do not include among 
the nuisance parameters the fine structure constant. This is because here we use the Fermi constant, mz and mw as 
input parameters, yielding a cm as output. 

Using notation consistent with previous analyzes we define our eight dimensional basis parameter set as 

m=((9,V) (11) 

which we will be scanning simultaneously over. For each choice of m a number of colliders or cosmological observables 
are calculated. These derived variables are denoted by £ = ^2> • • •)> which are then compared with the relevant 
measured data d. 

The quantity we are interested in is the posterior probability density function, (or simply posterior) p(m\d) which 
gives the probability of the parameters after the constraints coming from the data have been applied. The posterior 
follows from Bayes' theorem, 

P(m\d) = (12) 

where p(d\^), taken as a function of £ for fixed data d, is called the likelihood (where the dependence of £(m) is 
understood). The likelihood is the quantity that compares the data with the derived observables. 7r(m) is the prior 
which encodes our state of knowledge of the parameters before comparison with the data. This state of knowledge 
is then updated by the likelihood to give us the posterior. p{d) is called the evidence or model likelihood, and in 
our analysis can be treated as a normalization factor and hence is ignored subsequently for an example of how the 
evidence can be used for model comparison purposes. 

As our main prior we take very wide ranges of the CNMSSM parameters as given in Table HI although we have 
performed a number of additional scans which will be discussed below. We adopt a flat prior in logroi/2, log too > 
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CNMSSM parameters 
50 < mx/2 < 4 TeV 
50 < mo < 4 TeV 
\Ao\ < 7 TeV 
2 < tan/3 < 65 
10~ a < A < 0.7 
SM (nuisance) parameters tp 
160 < M t < 190 GeV 
4 < mt,(m b ) W < 5 GeV 
0.10 < q s (M z ) w < 0.13 

TABLE I: Initial ranges for our basis parameters m = (6,ip). 



tan j3 and A. Following Ref. [20| . we call this choice a log prior, as opposed to a completely /?ai prior used in some of 
our earlier analyzes where all the basis parameters are scanned with a flat measure. 

As before [20, our rationale for this choice of priors is that they are distinctively different. One reason why 
we apply different priors to soft mass parameters only is that they play a dominant role in the determination of the 
masses of the superpartners and Higgs bosons. Another important reason is that flat priors suffer from the "volume 
effect" by putting effectively too much emphasis on larger values of scanned parameters, while the log prior is more 
suitable for exploring smaller values of both m^ 2 and Too which are anyway more natural in effective low-energy 
SUSY models. Therefore the choice of log priors appears to be actually more suitable for revealing the structure of 
the model's parameter space, similarly as in the CMSSM [5l| and the NUHM (20l |. 

For the nuisance parameters we use flat priors (although this is not important as they are directly constrained by 
measurements) and apply Gaussian likelihoods representing the experimental observations (see table HI]), as before [lg, 

El MM- 

We compute our mass spectra and observable quantities using the publicly available NMSSMTools (version 2.1.1) 
that includes NMSPEC with a link to Micromegas; for details see Ref. |26[. We list the observables that the current 
version of NMSPEC, as linked with statistical subroutines available in SuperBayeS allows us to include in the likelihood 
function in Table [TTTT The relic density fl x h 2 of the lightest neutralino is computed with the help of Micromegas, 
which is also linked to NMSPEC. We further use the same code to compute the cross section for direct detection of 
dark matter via its elastic scatterings with targets in underground detectors but do not include it in the likelihood 
due to large astrophysical uncertainties. 

The likelihoods for the measured observables are taken as Gaussian with mean fi, experimental and theoretical 
errors (see the detailed explanation in Refs. 0, 53] )• In the case where there only an experimental limit is available, 
this is given, along with the theoretical error. The smearin g ou t of bounds and combination of experimental and 
theoretical errors is handled in an identical manner to Refs. [16l. [i~7j . with the notable exception of the Higgs mass 
and LEP limits on sparticle masses, which are calculated as a step function with values of the cross section times 
branching ratio (in the case of the Higgs) or mass that are within two standard deviations of the experimental limit 
being accepted. Finally, any points that fail to provide radiative EWSB, give us tachyons or the LSP other than the 
neutralino are rejected. 

As our scanning technique we adopt a "nested sampling" method [H} as implemented in the MultiNest 53] algo- 
rithm, which computes the Bayesian evidence primarily but produces posterior pdfs in the process. MultiNest provides 
an extremely efficient sampler even for likelihood functions defined over a parameter space of large dimensionality 
with a very complex structure. (See, e.g., Refs. (20l. [5l|.) This aspect is very important for the model analyzed here 
since the 8-dimensional likelihood hyperspace is fragmented and features many finely tuned regions that are difficult 
to explore with conventional fixed grid, random scan or even MCMC methods. For a comparison of CMSSM posterior 
maps obtained with a Metropolis-Hastings MCMC algorithm [l6|, [l?], EH and the MultiNest algorithm see Ref. [fOJ . 

As we are using nested sampling in this study, the issue of stopping criteria is handled differently from the MCMC 
case used in some earlier paper [lq . ITU [Lsj . Our treatments follows closely that presented in Appendix A of Ref. [5l[ . 
In nested sampling one is calculating the Bayesian evidence, defined by, 

Z = p(d) = f C{X)dX, (13) 
Jo 

where C is the likelihood and X the prior volume. One can get the posterior in a nested sampling scan, but the 
principle value calculated is the evidence. The stopping criteria takes into account that in general one is proceeding 
through shells of increasing iso-likelihood contours, with the set of "live points" drawn from within these contours. 
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SM (nuisance) parameter 


Mean value Uncertainty 
/! a (cxper.) 


Ref. 


M t 

m b (m b ) MS 
a B (M z ) MS 


172.6 GeV 1.4 GeV 
4.20 GeV 0.07 GeV 
0.1176 0.002 


f34] 
[35] 
[35] 



TABLE II: Experimental mean /j, and standard deviation a adopted for the likelihood function for SM (nuisance) parameters, 
assumed to be described by a Gaussian distribution. 



Observable 


Mean value 




Uncertainties 


Ref. 








a (exper.) 


t (theor.) 




8(g-2)„ x 10 10 


29.5 




8.8 


1 


[37] 


BR(B -> X s ~i) x 10 4 


3.55 




0.26 


0.21 


[37] 


BR(B U -> tv) x 10 4 


1.32 




0.49 


0.38 


[38] 


tl x h 2 


0.1099 




0.0062 


0.1{l x h 2 


[39] 




Limit (95% CL) 






t (theor.) 


Ref. 


BR(B S -» 


< 5.8 x 10" 8 






14% 


[40] 




As implemented in 


NMSSMTool. 






[26] 


sparticle masses 


As implemented in 


NMSSMTool. 






[26] 



TABLE III: Summary of the observables used in the analysis. Upper part: Observables for which a positive measurement has 
been made. 5(g — 2) tl denotes the discrepancy between the experimental value and the SM prediction of the anomalous magnetic 
moment of the muon (g — 2) M . For central values of the SM input parameters used here, the SM value of BR(B —* Xsj) is 
3.11 x 10 -4 , while the theoretical error of 0.21 x 10 -4 includes uncertainties other than the parametric dependence on the 
SM nuisance parameters, especially on M t and a a (Mz) MS ■ For each quantity we use a likelihood function with mean n and 
standard deviation s — \/<r 2 + r 2 , where a is the experimental uncertainty and r represents our estimate of the theoretical 
uncertainty (see Ref. [l(| for details). Lower part: Observables for which only limits currently exist. The likelihood function is 
given in Ref. [r|, including in particular a smearing out of experimental errors and limits to include an appropriate theoretical 
uncertainty in the observables in BR(B S — > fi + jj,~). The limit on the light Higgs mass mj, is applied in a simplified way, see 
text for details. 

One can then define the stopping criterion as taking the maximum likelihood point in the set of live points, see 
Ref. [52| (£ ma x) and calculating the maximum change to the evidence it could make, 5Zi = £, max Xi. Once this value 
goes below a specified value (we take SZ < 0.5) the run is terminated. 

IV. PROBABILITY MAPS OF CNMSSM PARAMETERS AND OBSERVABLES 

In this section we present our numerical results from global scans of the CNMSSM parameter space. We begin 
with the CNMSSM parameters and next show probability maps for several observables, including, in turn, the Higgs 
bosons, some superpartners and other collider signatures, and finally dark matter cross sections. 

To begin with, in Fig. [T] we plot joint 2D relative probability density functions (pdfs) for some combinations of the 
CNMSSM parameters in our default case as given in Table Q] and taking the log prior, as explained above. In this, 
and figures below showing 2D pdfs, the inner (outer) contours delineate the 68% (95%) total probability regions and 
the color code is given in the bar at the bottom. 

First, we can see that higher probability regions for all the parameters but mo are confined well within the assumed 
priors and show clear high probability peaks. Focusing on the left panel in the plane spanned by m 1( / 2 and mo, we can 
see some prominent features: a rather strong preference for the stau coannihilation region of m 1 / 2 > mo ~ 0.5 TeV, 
although the 68% total probability region extends to larger m because of the pseudoscalar funnel effect contribution 
to Q x h 2 , and to much larger values of the parameter of the focus point (FP) region [42],[43j]. Green triangles indicate 
some of the best fit points. The wedge of m 1 / 2 3> mo is disallowed because of charged LSP (normally the stau). 

Examining the other panels of Fig. Q] one can see a rather strong preference for large tan (3 (as in the CMSSM) , 
although small values below some 15 are also favored. Aq appears to have two distinct branches of opposite sign, with 
Aq = 0, although not excluded, proving to be hard to find solutions for, with some resemblance to the CMSSM. In 
some sense the latter may even be as favored as large positive values, as indicated by one of the best fit points. 

Clearly, these are familiar features of the CMSSM, as one can see by comparing the high probability regions of the 
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FIG. 1: The 2D relative probability density functions in the planes spanned by the CNMSSM parameters m^, mo, tan/9 and 
Aq for the log prior. The pdfs are normalized to unity at their peak. The inner (outer) blue solid contours delimit regions 
encompassing 68% and 95% of the total probability, respectively. All other basis parameters, both CNMSSM and SM ones, in 
each plane have been marginalized over (i.e., integrated out). Blue dots denote some best fit points. 
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CNMSSM in Fig. [T] with analogous figures for the common set of parameters shared with the CMSSM, as shown in 
Fig. 13 of Ref. [lyp (obtained with the NS scan), or with Fig. 1 of Ref. [H[ (obtained with the MCMC scan). 

The similarity of the high probability regions of m.1/2 and mo in both models suggests the parameters of the 
CNMSSM tend to favor the decoupling limit, A — > 0. In Fig. [H we show 2D pdfs of A with the CMSSM-like 
parameters. One can see that in general A prefers to be small, which leads to a statistical preference for a very 
CMSSM-like behavior. Large values of A, bigger than around 0.6 are disfavored due to a Landau pole in the running 
of A. At "intermediate" values, 0.1 < A < 0.6, the constraints become weaker but there remain problems with 
tachyons, seen most clearly in the (A, k) plane, shown in Fig. [31 which also shows how both parameters are rather 
closely correlated, k oc A, and favor small values, towards the decoupling limit. The region with k 3> A is disfavored 
by the presence of tachyonic CP-odd scalars, and similarly with CP-even scalars for k <C A. The preference for low A 
could be due to the fact that there are fewer tachyonic directions in the potential close to the decoupling limit A — > 
(for this to happen A < 0.1 is sufficient). For a more detailed discussion, see Ref. [28j |. 

In conclusion, the parameters of the CNMSSM seem to favor the CMSSM limit since at small values of A < 0.1 
it is simply much easier to find physical solutions. This has been confirmed with exploratory runs with only the 
requirement of correct EWSB and a neutralino LSP enforced. Indeed, a flat (and nor a log) prior in A was chosen 
so as not to emphasize low values of the parameter and instead to "force" the scan away from the decoupling limit. 
Even with the flat prior, however, the preference for small A remains strong. For A below 0.1 there are more solutions 
because for instance tachyons are less of a problem (due to less mixing with singlets) and also we are further away 
from the Landau pole region. 

It is also instructive to show 2D probability maps of the effective \x parameter vs. some of the CNMSSM parameters. 
This is presented in Fig. [5] We can see that the consistency of the model and the applied set of constraints favor /i 
below some 2 TeV, the range comparable to Msusy, as expected. In other words, in the CNMSSM the \i problem 
is solved without any need for additional fine tuning of parameters. We can also see an interesting correlation with 
m 1 / 2 but not with the other C(N)MSSM parameters. This is caused primarily by the CP-odd Higgs a\ funnel effect 
and the fact that its mass is correlated with fi. These features are presented in Fig. El Finally, in Fig. [6] we present 
ID pdfs of several key parameters which show more clearly their high probability ranges. 

Prior dependence is often an issue in Bayesian statistics and needs to be addressed. Even the CMSSM tends to be 
under-constrained which leads to a fairly strong prior dependence [5l| . and in the NUHM, with two extra parameters, 
the situation becomes worse (20j . In the left panel of Fig. [7] we show a 2D pdf in the (?tt-i/2,™o) plane assuming 
the flat prior in all the CNMSSM parameters. By comparing with the analogous panel of Fig. [1] we indeed see a 
substantial shift in the high probability region to larger values, as typical for the flat prior due to the volume effect. 

A related issues is that of the assumed range of input parameters. We have already seen in the left panel of Fig. [7] 
that mo was not well confined to the assumed range below 4 TeV. In order to examine this, in the right panel 
of Fig. [7| we show a 2D pdf in the (m 1 / 2 ,m ) plane with greatly extended ranges of both parameters (50 GeV < 
rai/2,mii < 10 TeV) and taking the log prior. Clearly, there is basically no cap on the 95% total probability range 
in both parameters, although of course such large values of soft mass parameters can hardly be considered as well 
motivated in effective low-energy SUSY models. 

We have already emphasized that the high probability regions of the crucial parameters m 1 / 2 and mo in the 
CNMSSM are quite similar to the well studied CMSSM case. One of the key features of the CMSSM is that the 
neutralino LSP is mostly a bino, except for the FP region of large mo, where a larger admixture of the higgsino 
is present. In the CNMSSM, with an additional singlet field whose mass is controlled by ks, the picture could in 
principle be very different from the CMSSM. We examine this in Fig. [5J where we separately show the regions where 
the LSP is mostly gaugino (Z g = Z\ x + Z\ 2 > 0.7) (left panel), doublet higgsino (Z h = Zf 3 + Zf 4 > 0.7) (middle left 
panel), the mixed region (0.3 < Z g ,Zh < 0.7 and Z s = Z X5 < 0.5) (middle right panel), as well as mostly singlino 
(Z s > 0.5) (right panel). (The wino component is always negligible and below we will show only the bino fraction 
Zi, = Zfi-) We can see that in an overwhelming fraction of cases the LSP still remains predominantly bino-like. This 
is in agreement with the panel showing ks in Fig. [6] where small values of the product tend to be strongly disfavored. 
We expose those non-CMSSM like cases in Fig. Clearly, points corresponding to singlino LSP cases exist but are 
rare. 

Similarly to the CMSSM, the key observables shaping the high probability regions of the model are: fl x h 2 , BR(B — ► 
Xsl)i ${g ~ 2) M and the light Higgs mass mj. Their ID pdfs are presented in Fig. [TDJ We can see that, apart from 
5{g— 2)u, they reproduce the likelihood rather well, especially for our default log prior (long-dashed red curve), although 
even for the flat prior the fit is not bad, except for S(g — 2)„, again with much resemblance to the CMSSM (5l| . 

The sharp cutoff in m/, at the LEP limit results from the approximation that we have adopted, as mentioned 
above. Given the complexities of the Higgs sector and a much larger list of possible decay channels, we have taken 
an approximate approach of setting the likelihood to one or zero for points that are accepted or excluded by LEP 
data using the NMSPEC code. A more correct would would have allowed for some "tail" at lower masses, as in the 



CMSSM |17|, and would have expanded allowed ranges of m l / 2 and mo towards smaller values but this would have a 
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FIG. 2: The 2D relative probability density functions in the planes spanned by A and the other CNMSSM parameters for the 
log prior. The pdfs are normalized to unity at their peak. The inner (outer) blue solid contours delimit regions encompassing 
68% and 95% of the total probability, respectively. All other basis parameters, both CNMSSM and SM ones, in each plane 
have been marginalized over. Blue dots denote some best fit points. 



limited effect on this analysis in which we are mostly interested in presenting mainly global features of the CNMSSM. 
Nevertheless, it is worth stressing that in an exploratory run with the LEP limit smeared out we have found only a 
few points for which the singlet component is barely enough to escape the limit of 114.4 GeV. In this sense we expect 
that a full analysis would have basically reproduced the case of the CMSSM [TtJ ■ 

The Higgs sector of the NMSSM contains three CP-even bosons hi = h, hi and h%, two CP-odd ones a\ and <Z2, 
as well as a pair of the charged Higgs H ± . Fig. [TT1 shows the relative ID pdfs of hi (left panel), a% (middle panel) 
and H^ 1 (right panel) for both our default log prior (long-dashed red) and the flat prior (dotted blue) for comparison. 
We can see that the prior dependence is not very strong, with the fiat prior favoring larger values, as usual. Also, by 
comparing with Fig. 5 of Ref. 17j, we can see that the pdfs are quite similar to the analogous states in the CMSSM. 

In Fig. [12] we present the relative ID pdfs for several superpartners, in a fashion similar to the previous Figure. 
Again, we can see that the log prior gives somewhat lower ranges of masses, especially for the scalars, which primarily 
depend on mo and that the distributions are rather similar to the corresponding ones in the CMSSM; compare, eg, 
Fig. 17 of Ref. [IB]. It is clear that there will be a rather mixed chance of detecting those states at the LHC. For 
example, with the gluino to be probed up to some 2.7 TeV, nearly the whole range will be tested with the log prior, 
but much less so with the flat prior. The scalars, on the other hand, will be much more challenging for both priors. 
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FIG. 3: The 2D relative pdfs in the plane of (A, «) for the log prior. 



Parameter 


Best fit (log) 


Best fit (flat) 


mi/2 


fOf GeV 


478 GeV 


mo 


404 GeV 


632 GeV 


A 


-165 GeV 


1.20 TeV 


tan/3 


12.9 


42.4 


A 


0.009 


0.0252 


A 4 


547 GeV 


672 GeV 


m 01 


274 GeV 


476 GeV 




0.093 


0.094 


br(b -> X s7 ) 


3.10 x 10~ 4 


3.27 x 10" 4 


Bi?(B s -» /Lt+jT) 


2.8 x 10" 9 


1.6 x 10~ 8 


BR{B U -> tv) 


1.28 x 10 -4 


0.93 x 10" 4 


«(ff-2)„ 


16.9 x 10" 10 


14.8 x 10" 10 


m h 


114.4 GeV 


114.3 GeV 


m x 


164 GeV 


263 GeV 


m ± 


309 GeV 


491 GeV 




950 GeV 


1.45 TeV 


X2 


9.6965 


9.4635 



TABLE IV: A table showing the values of various parameters for the best fitting point in both the log and flat prior case. 



In Tabic [TVl we list the best fit values for a number of observables for our default log prior and also, for comparison, 
for the flat prior. We stress, however, that the log prior appears more appropriate for exploring unified low-energy 
SUSY models, as already emphasized above. 

Finally, we move to discussing the model's predictions for the detection of the lightest neutralino assumed to be 
the DM in the Universe in direct detection searches via its elastic scatterings with targets in underground detectors. 
We follow the same procedure and formalism as previously in pH Il8l . [5l| . The underlying formalism can be found in 
several sources, e.g., in Refs. [54l. l55l. l56l . 

In Fig. [13] we present the 2D posterior pdfs in the usual plane spanned by the spin- independent cross section 1 
and the neutralino mass m x . The left (right) panel corresponds to the log (flat) prior. For comparison, some of the 
most stringent 90% CL experimental upper limits are also marked 0, |45|, |4(| |48|, Ei| , although they have not been 
imposed in the likelihood, as before in our studies of the CMSSM, because of substantial astrophysical uncertainties, 
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FIG. 4: The 2D relative probability density functions in the planes spanned by fi and the CNMSSM parameters that are the 
same for the CMSSM for the log prior. The pdfs are normalized to unity at their peak. The inner (outer) blue solid contours 
delimit regions encompassing 68% and 95% of the total probability, respectively. All other basis parameters, both CNMSSM 
and SM ones, in each plane have been marginalized over. Blue dots denote some best fit points. 



especially in the figure for the local density. 

Several key features can be seen in Fig. [131 Firstly, the prior dependence is not very strong for both 68% and 
95% total probability regions, which is encouraging. It does not affect much the banana-shape high-probability region 
which corresponds to the Higgs funnel and the stau coannihilation regions. The horizontal branch of 1 ~ 7 x 10~ 8 pb 
is more affected because it corresponds to the focus point region of large mo. Next, the overall shape rather closely 
resembles the case of the CMSSM, see, e.g., Fig. 18 of Ref. [Hj or Fig. 13 of Ref. 16]. (The slight upwards shift 
in <7p J results from changing the code from DarkSusy to Micromegas.) It does, on the other hand, differ from the 
predictions of the NUHM which features an additional higgsino-like region at m x ~ 1 TeV; see Fig. 12 of Ref. [20]. 

Independently of the prior, basically the whole 68% and 95% total probability regions are likely to be within the 
planned reach of 10~ 10 ± of future 1-tonne detectors. Some of the currently operating detectors are already probing 
some of the high probability regions, and with a "modest" improvement down to ~ 10~ 8 pb, they will be testing some 
of the most likely cross sections. 
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FIG. 5: The 2D relative probability density functions in the plane of (m x ,m ai ) (left panel) and (p,m ai ) (right panel) for the 
log prior. The pdfs are normalized to unity at their peak. The inner (outer) blue solid contours delimit regions encompassing 
68% and 95% of the total probability, respectively. All other basis parameters, both CNMSSM and SM ones, in each plane 
have been marginalized over. 
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FIG. 6: The ID relative probability densities for all the CNMSSM parameters, plus k, fi, As and ks. 



V. CONCLUSIONS AND SUMMARY 



The Next-to-Minimal Supersymmetric Standard Model solves the ^-problem of the MSSM but, without grand 
unification, both models suffer from a large number of parameters. The constrained versions of both models are in 
this respect much more well-motivated. Because of the additional singlet superfield present in the CNMSSM, the 
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FIG. 7: The same as in fig. \T\ but for a flat prior (left panel) and a log prior with a greatly extended range, 50 GeV < 
m i/2i m o < 10 TeV (right panel). 



resulting phenomenology in the Higgs and neutralino sectors is considerably richer. Therefore, a prior one could 
expect that the models could be distinguished in experimental tests. 

The global exploration of wide ranges of CNMSSM parameters and a Bayesian analysis show that, from the 
statistical point of view, this is not the case. The coupling A strongly favors as small values as possible, in other 
words it tends towards the decoupling regime in which one recovers the CMSSM plus the basically decoupled singlet 
Higgs and the singlino. As a result, Higgs and superpartner mass spectra also tend to resemble those of the CMSSM, 
as does the cross section for direct detection of neutralino dark matter. Nevertheless, we have identified a limited 
number of cases where the LSP is indeed singlino-dominated, but statistically they are not very significant. 

In conclusion, the CNMSSM is, for the most part, testable at the LHC and in dark matter searches, which is 
certainly encouraging. On the other hand, should a CMSSM-like signal be detected, it is likely to be very challenging 
to distinguish between the two models. 
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m ± (upper middle panel), the gluino rrig (upper right panel), the lighter stop t\ (lower left panel), left squark <jl (lower middle 

panel) and the lighter stau t\ (lower right panel). In each panel we show the posterior for the flat prior (dotted blue) and the 
log prior (long-dashed red). 
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